      SUBROUTINE SETUP
      RETURN
      END
C   HS32
      BLOCK DATA
      INCLUDE 'BLOCO.INC'
      INTEGER I,J
      DATA NAME/'HS32'/
      DATA X(1),X(2),X(3)/ .1D0,.7D0,.2D0/
      DATA N/3/ , NH/1/ , NG/ 4/
      DATA DEL0/1.00D0/ ,TAU0/1.D0/ ,TAU/.1D0/
      DATA (GUNIT(1,J),J=2,4)/3*1/,(GUNIT(3,I),I=2,4)/3*1/,
     F     (GUNIT(2,I),I=2,4)/1,2,3/
     F     ,((GUNIT(I,J),I=1,3),J=0,1)/-1,0,0,-1,0,0/
     F     ,((GUNIT(I,J),I=1,3),J=5,5)/-1,0,0/
      END
      SUBROUTINE EF(X,FX)
      INCLUDE 'FUCO.INC'
      DOUBLE PRECISION X(NX),FX
      ICF=ICF+1
      FX=(X(1)+3.D0*X(2)+ X(3))**2+4.D0*(X(1)-X(2))**2
      RETURN
      END
      SUBROUTINE EGRADF(X,GRADF)
      INCLUDE 'FUCO.INC'
      DOUBLE PRECISION X(NX),GRADF(NX),XKLA1,XKLA2
      ICGF=ICGF+1
      XKLA1=X(1)+3.D0*X(2)+X(3)
      XKLA2=X(1)-X(2)
      GRADF(1)= 2.0D0* XKLA1+8.D0*XKLA2
      GRADF(2)= 6.D0* XKLA1-8.D0* XKLA2
      GRADF(3)= 2.D0* XKLA1
      RETURN
      END
      SUBROUTINE EH(I,X,HXI)
      INCLUDE 'FUCO.INC'
      INTEGER I
      DOUBLE PRECISION X(NX),HXI
      CRES(I)=CRES(I)+1
      HXI= 1.D0-X(1)-X(2)-X(3)
      RETURN
      END
      SUBROUTINE EGRADH(I,X,GRADHI)
      INCLUDE 'FUCO.INC'
      INTEGER I,J
      DOUBLE PRECISION X(NX),GRADHI(NX)
      CGRES(I)=CGRES(I)+1
      DO      100      J=1,3
      GRADHI(J)=-1.D0
 100  CONTINUE
      RETURN
      END
      SUBROUTINE EG(I,X,GXI)
      INCLUDE 'FUCO.INC'
      INTEGER I
      DOUBLE PRECISION X(NX),GXI
      IF ( GUNIT(1,I+NH) .EQ. -1 ) CRES(I+NH)=CRES(I+NH)+1
      IF(I .GT. 3)       GOTO 100
      GXI=X(I)
      RETURN
  100 CONTINUE
      GXI=6.D0*X(2)+4.D0*X(3)-X(1)**3-3.D0
      END
      SUBROUTINE EGRADG(I,X,GRADGI)
      INCLUDE 'FUCO.INC'
      INTEGER I,J
      DOUBLE PRECISION X(NX) ,GRADGI(NX)
      IF(I .GT. 3)       GOTO 100
      DO      10      J=1,3
      GRADGI(J)=0.D0
  10  CONTINUE
      GRADGI(I)=1.D0
      RETURN
  100 CONTINUE
      GRADGI(1)=-3.D0*X(1)**2
      GRADGI(2)= 6.D0
      GRADGI(3)= 4.D0
      RETURN
      END


